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Abstract 

According  to  the  Mars  Express  mission,  the  MARSIS 
primary  scientific  objectives  are  to  map  the  distribution  of 
water,  both  liquid  and  solid,  in  the  upper  portions  of  the 
crust  of  Mars.  Three  secondary  objectives  are  also  defined 
subsurface  geologic  probing,  surface  characterization,  and 
ionosphere  sounding. 

In  order  to  obtain  the  primary  objectives  the  Radar  Sounder 
design  was  based  on  the  Ice/water  interface  and  Dry/ice 
interface  scenario:  defining  the  material  composition  of  the 
first  layers  and  porosity  and  the  pore  filling  materials. 
Concerning  the  surface,  we  have  characterized  the 
geometric  structure  in  terms  of  a  large-scale  morphology, 
on  which  a  small-scale  geometric  structure,  due  to  rocks,  is 
superimposed,  taking  into  account  also  that  recently  the 
structure  of  the  planets  surface  was  described  by  means  of 
fractals  and  in  particular  the  new  MARS  surface  models 
obtained  by  processing  of  the  MOLA  data. 

According  to  these  models,  this  paper  provides  a  description 
of  the  operational  planning  approach  and  expected 
performances  of  MARSIS. 

1  INTRODUCTION 

MARSIS  is  a  low-frequency  nadir-looking  pulse  limited 
radar  sounder  and  altimeter  with  ground  penetration 
capabilities,  this  radar  can  effectively  operate  at  any  altitude 
lower  than  800  km.  In  order  to  maximize  the  penetration 
capabilities  of  the  transmitted  pulse  MARSIS  must  operate 
at  a  frequency  as  low  as  possible  as  carrier  frequencies 
taking  into  account  the  expected  values  of  the  plasma 
frequency  in  the  Martian  ionosphere:  1.8,  3,  4  and  5  MHz 
was  selected.  Moreover  the  requirement  of  fine  range 
resolution  entails  a  relatively  large  transmitted  bandwidth  (1 
MHz),  so  that  MARSIS  will  operate  with  a  very  high 
fractional  bandwidth:  a  1  MHz  bandwidth  allows  a  vertical 
resolution  of  150  m  in  vacuum  which  corresponds  to  50-100 
m  in  the  subsurface,  depending  on  the  e.m.  wave 
propagation  speed  in  the  crust. 

According  to  the  well  known  principle  of  operation  of  a 
subsurface  radar,  a  short  pulse  of  e.m.  energy  transmitted  by 
the  antenna  impinging  on  the  top  of  the  Mars  surface 
produce  a  first  reflection  echo,  which  propagates  backward 
to  the  radar,  moreover,  thanks  to  the  long  wavelengths 
employed,  a  significant  fraction  of  the  e.m.  energy 
impinging  on  the  surface  is  transmitted  into  the  crust  and 
propagates  downward.  Additional  reflections,  due  to  the 
subsurface  dielectric  discontinuities  would  occur  and  the 


relevant  echoes  would  propagate  backward  through  the  first 
layer  medium  and  then  to  the  radar  generating  further  echo 
signals,  much  weaker  than  the  front  surface  signal.  The  time 
delay  of  the  echo  can  be  converted  to  a  depth,  assuming  that 
the  propagation  speed  of  the  medium  is  known,  and  the 
intensity  of  the  reflection  can  be  analyzed  to  estimate  the 
reflectivity  at  the  interface  and  the  attenuation  properties  of 
the  intervening  layers. 

In  order  to  optimize  the  radar  design,  according  to  the 
suggestions  from  planetary  scientists,  we  have  selected 
some  models  of  subsurface  porosity  profile  and 
composition,  as  described  in  the  following.  Moreover  the 
two  most  likely  scenarios  representing  the  relevant 
interfaces  are  the  following: 

Ice/water (I/W)  interface  -  in  this  scenario  the  pores  are 
filled  with  ice  from  the  surface  down  to  a  depth  below 
which  liquid  water  is  stable  and  becomes  the  pore-filling 
material.  The  change  of  the  pore-filling  material  causes  a 
discontinuity  of  the  overall  dielectric  constant,  which  can  be 
detected  by  the  radar  sounder. 

Dry/ice(D/I)  interface  -  here  the  pore-filling  material  is 
considered  to  be  gas  or  some  other  vacuum-equivalent 
material  up  to  a  depth,  below  which  ice  fills  the  pores. 
Hence  the  interface  to  detect  is  between  dry  regolith  and 
ice-filled  regolith. 

We  can  notice  that  a  dynamic  range  of  50-60  dB  can  be 
managed  by  the  instrument,  but  there  are  several  factors  that 
can  strongly  reduce  the  detection  dynamic  mainly  the  noise 
and  the  surface  clutter  In  fact  due  to  the  orbital  sounding 
geometry,  the  off  nadir  surface  return  will  be  received  at  the 
same  time  of  subsurface  echoes.  Then,  in  case  of  rough 
surface,  the  surface  echoes  may  happen  to  hide  weak 
subsurface  features  and  the  detection  of  discontinuities  in 
the  crust  could  be  indered  by  surface  clutter.  As  more  the 
surface  is  smooth  as  more  rapidly  the  returns  from  the 
surface  at  off-nadir  angles  will  drop  off. 

Hence  it  is  important  to  evaluate  the  penetration  depth  that 
is  possible  to  reach  according  to  the  transmitted  wavelength 
the  dielectric  characteristics  of  the  crust  and  the  surface 
scattering  behavior  as  a  function  of  the  surface  topography, 
which  uses  synthetic  aperture  techniques  and  a  secondary 
receiving  antenna  to  isolate  subsurface  reflections. 

Moreover  different  techniques  are  envisaged  to  increase  the 
detection  performance  against  surface  clutter.  First  of  all  the 
Doppler  Azimuth  Processing  significantly  reduces  the 
surface  echoes  coming  from  along  track  off  nadir 
reflections:  the  improvement  on  the  overall  surface  clutter 
attenuation  should  be  of  the  order  of  10  dB. 
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Then  the  instrument  secondary  antenna ,  a  monopole 
oriented  along  the  nadir  axis,  will  receive  mostly  the  off- 
nadir  surface  returns,  that  could  be  thus  subtracted  by  the 
primary  antenna  composite  signal,  further  reducing  the 
surface  clutter  level  (about  ~15dB). 

Finally  echo  profiles  collected  at  two  different  frequencies 
can  be  processed  to  separate  the  subsurface  reflections, 
which  are  strongly  dependent  on  the  frequency,  from  the 
surface  reflections,  which  are  mostly  frequency  independent 
(the  achieved  improvement  can  reach  10-15  dB)  and  to  have 
information  of  the  first  layer  attenuation. 


2  SURFACE  AND  SUBSURFACE  SCATTERING 
MODELS 

In  order  to  estimate  the  reflectivity  at  the  interface  and  the 
attenuation  properties  of  the  intervening  layers  and  to 
assess  the  interface  detection  performance  it  is  required  to 
evaluate  the  back  scattering  cross  sections  of  concurrent 
echoes  coming  from  the  surface  and  subsurface  layers: 
Os  =  r  •f,(H1>s1,x)  oss  =  r,  r„(ii„.s„./v) 
being  Ts  and  TSs  the  Fresnel  reflectivity  terms,  which  deal 
with  the  surface  and  subsurface  dielectric  properties  and  fs 
and  fss  the  scattering  terms,  which  deal  with  the  geometric 
structure  of  the  surface  and  subsurface. 

2. 1  Fresnel  coefficients  and  subsurface  attenuation 
In  Tab.l  are  listed  the  values  of  dielectric  constant  of  the 
materials  which  can  be  considered  as  end  members  of  the 
range  in  which  the  first  layer  of  the  Martian  surface 
materials  may  vary.  The  dielectric  properties  of  the  water 
and  ice  filling  the  pores  are  listed  in  the  Tab. 2. 

It  should  be  noted  that,  since  porosity  depends  on  the  depth, 
so  will  do  the  effective  dielectric  constants  of  the  mixture. 
In  order  to  evaluate  the  mixture  dielectric  constants,  we 
have  used  the  Maxwell-Gamett  model 
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In  (2.1.2)  K  is  the  decay  constant  that,  for  Mars,  can  be 
computed  by  scaling  the  measured  Lunar  decay  constant  for 
the  ratio  between  the  Lunar  and  Martian  surface 
gravitational  accelerations,  under  the  assumption  of 
comparable  crust  densities:  the  resulting  value  for  Mars  is 
K=2.8  km. 


Using  the  eq.(2.1.1),  the  effective  dielectric  constant  of  the 
mixtures  which  have  been  selected  to  represent  the  Martian 
porous  regolith  can  be  evaluated,  as  a  function  of  the  depth, 
and  of  the  surface  porosity.  Therefore  the  well  known 
Surface  Fresnel  reflectivity  for  nadir  incidence  (from  free 
space  propagation)  can  be  evaluated 
(2.1.3) 


rs  = 


2 
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being  erl(0)  the  real  dielectric  constant  of  the  crust 
evaluated  at  the  surface  (z=0). 


The  selected  models  make  the  surface  reflectivity  ranging 
between  -6.5  dB  and  -9.5  dB  for  the  I/W  scenario  and  -7  dB 
and  -12  dB  for  the  D/I  scenario. 
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20% 

I 

-9.5 

-9 

-12 

-9.5 

II 

-7.5 

-6.5 

-9 

-7 

III 

-8.5 

-7.5 

-10 

-8 

Tab. 3 


The  Fresnel  reflectivity  for  a  subsurface  layer  located  at  a 
depth  z’  can  be  expressed  as  follows 

z1 

-0.1  ja{z)d2 

(2.1.4)  ^,=44-^,)  10  • 

being  R2i2,z’  the  reflection  coefficient  of  such  interface: 

(2-!-5)  <z,= 

By  the  defined  crust  materials,  we  have  obtained  for  the  two 
scenario: 

IAV  ->  R2i2,z-=  (-10-5—18)  -  3dB/Km 

(2.1.6) 

D/I  R2i2,z’=  (-20+-28)  -  3dB/Km 
where  the  first  term  is  related  to  the  first  layer  porosity  and 
the  last  one  is  due  to  the  decreasing  of  the  porosity  with  the 
depth.  We  can  notice  that  the  R2i2,z’  is  near  crust  material 
independent. 

Concerning  the  two  way  attenuation: 

(2.1.7)  adB/m{z)  =  1.8-1  O'7 f0 FF) tan 8{z) 

we  have  obtained,  as  a  good  approximation,  the  functions 
detailed  in  Tab.4. 
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(0.8+0.  lz) 

(1.3+0.05z) 

(0.9+0.  lz) 

(1.3+0.05z) 

ii-l/m 

(3.7+0.54*z) 

(6+0.25  *z) 

(4.3+0.54*z) 

(6.3+0.25*z) 

II-2 

(8.5+1.  l*z) 

(13.3+0.5*z) 

(9.5+1.  l*z) 

(13.7+0.5*z) 

Tab.4  (z  depth  in  Km) 


2.2  Surface  backs cattering  models 

At  the  beginning  of  our  analysis  we  modeled  the  Martian 
surface  back  scattering  considering  two  main  terms:  the 
first,  called  large  scale  scattering  contribution ,  resulting 
from  gentle  geometrical  undulations  of  the  surface  on  a 
scale  of  many  hundreds  to  thousands  meters,  whereas  the 
second,  called  small  scale  scattering  contribution ,  giving 
reason  of  the  fast,  slight  variations  of  the  surface  height  over 
an  horizontal  scale  of  some  tenths  of  meters.  A  plausible 
range  for  the  parameters  describing  the  surface  geometry  is 
listed  in  the  Tab. 5. 


LARGE  SCALE  MODEL  | 

SMALL  SCALE  MODEL 

RMS  slope  (ms) 

Correlation  Length 

RMS  Slope 

RMS  height 

0.01  -0.1  rad. 
(0.57°- 5. 7°) 

200  -  30000  m 

0.1  -  0.6  rad. 
(5.7°  -34.3°) 

0.1  - 1  m 

Tab.5:  Summary  of  the  values  range  for  the  geometric 
parameters  of  the  surface 
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From  the  statistical  point  of  view,  the  surface  height  was  be 
modeled  as  a  Gaussian  random  process,  being  Oh  the  total 
surface  RMS  height.  At  a  second  order  level  we  suppose  (as 
usual)  an  isotropic  correlation  coefficient. 

Simple  approximate  methods  can  be  applied  for  surfaces 
which  present  a  unique  roughness  scale,  with  either  a  big 
correlation  length  (gently  undulating  surface),  or  a  very 
small  rms  height  (slightly  rough  surface)  compared  to  the 
incident  wavelength.  Specifically,  the  Kirchhoff  method  can 
be  applied  for  gently  undulating  surfaces,  which  respect  the 
tangent  plane  conditions  and  the  Small  Perturbation  Method 
can  be  applied  to  slightly  rough  surfaces.  The  classical 
studies  on  the  validity  conditions  of  these  two  models  [1] 
have  been  recently  updated  [2],  and  regions  of  validity 
currently  defined  are  reported  in  Fig.l.  It  is  evident  from 
these  figures  that  the  Kirchhoff  approximation  can  be  used 
to  evaluate  the  large  scale  backscattering  contribution 
whereas  the  Small  perturbation  method  can  be  used  for  the 
small  scale  contribution. 


Fig.l:  KIRCHHOFF  and  SMALL  PERTURBATION  approximation 
validity  conditions  {k=2nlX) 


2.2.1  Large  Scale  Contribution  evaluation 

Under  Kirchhoff  approximation  hypotheses  the  scattered 
electric  field  to  the  antenna  is  given  by  the  following 
integral  [1]: 

(2'2'1)  Es(t)  =L  j  fiaj) .  G(P)  •  Mp)  ■  •  em<'  ds-  e'md(0 

where  Ei=(l/R1)e"jkRi  is  the  incident  spherical  wave,  F(co)  is 
the  Fourier  transform  of  the  transmitted  pulsed  waveform,  A 
is  the  wavelength,  P  is  the  generic  point  on  the  surface, 
G(P)  is  the  antenna  normalized  gain  in  the  direction  of  the 
point  P  on  the  surface  S,  R(P)  is  the  Fresnel  coefficient, 
is  the  distance  from  radar  to  the  point  P,  Rj  is  the  unit 

vector  corresponding  to  the  path  from  the  radar  to  the  point 
P ,  n  is  the  unit  vector  normal  to  the  surface  in  the  point  P 
and  finally  k=2jt/X  is  the  wave  number. 

If  we  make  the  further  hypotheses  of  Isotropic  antenna 
pattern  and  surface  tilts  are  small  enough  to  allow  that  we 


confuse  the  local  normal  to  the  surface  with  the  vertical 
axis,  we  have  G(P)=1  and  n  •  Rx  ~  z  •  R{  =  COSO  ,  where 
0  is  the  incidence  angle  with  respect  to  the  vertical. 
Moreover  the  distance  from  the  radar  to  the  generic  point 
P(x,y,z)  can  be  approximated,  under  the  far-field 
approximation,  as: 

(2.2.2)  R  ~h-zcos  6  +  x2+y2 
1  2  h 

where  h  is  the  distance  from  the  Radar  to  the  mean  surface. 
With  all  these  assumptions  the  eq.  (2.2.1)  becomes: 

oo  •  /  x2+y2  2z-cos6?) 

<2-2-3 )  Es (t)  =  -2^6-  r(0) f (is  [ e' r  h'c+  c  > jccF(co)dco 
An -eh  {  J 


where  T —t - and  the  angle  of  incidence  has  been 

c 

supposed  to  undergo  no  significant  variation  within  the 
pulse  limited  integration  area,  and  has  been  therefore  taken 
out  of  the  surface  integral. 

The  average  scattered  power  at  time  T  can  be  evaluated 
taking  the  average  of  the  product  of  the  scattered  electric 
field  with  its  complex  conjugate: 

(eA)AW)  =  /  4r^v;  cos2  ojdcqjdaw  ■  c^F^Y ■ 
[An-chf 

(2-2-4)  ■  > 

Si  S2 

where  r,(0)  =  |/?(0)2  is  the  Fresnel  reflection  coefficient. 

The  average  on  the  height  z  can  be  expressed  as  a  function 
of  the  surface  two-dimensional  characteristic  function: 


/^(CD^-CQ^  )  -2oJiCOS2  6(of  +of-2p(-Y1  ,x2 ,  /! ,  72  ) 

(2.2.5)  <e  c  >=e 

The  solution  of  the  integral  is  not  easily  achievable  in  the 
general  case:  a  closed  form  has  been  derived  [4]  assuming  a 
gaussian  correlation  coefficient  (gaussian  roughness 
spectrum): 


(2.2.6)  p(£)=e  i2  <=> 


Pn 


which  corresponds  to  a  gaussian  shape  for  the  compressed 
pulse 

t2 

(2.2.7)  /(/)=  e  4cT?  cos  CO0t 


where  C00  is  the  carrier  frequency  and  aP  is  related  to  the 

0  ^7 

system  bandwidth  Bw :  G  . 

P  Bw 

These  hypotheses  lead  to  the  following  results  [4] : 

(2.2.8)  <4(rK(T)>=r,(0)A.Cos20  (pc  +  pna-pj 

where  Pc  (x)  is  the  coherent  (specular)  scattering 
component: 

^  (2  k  <7 /,]  cos  6  )2  T _ 

(2.2.9)  P  (t)=  — — e  ^  •  e  2ct^1+f) 

1  +  F 

while  Pnc  (%)  =  R nci (% ) ~Rnc2 (% )  is  the  non- coherent  (diffuse) 
scattering  component: 
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F  = 
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surface  parameter 

Therefore  the  cross  section  of  the  large  scale  surface  model 
is  given  by: 

(2.2.14) 


C7i(t)  = 


4^2  =r  (0M2  cos2  0  +  +Jp„cl  -P„c2) 


We  can  notice  that  the  maximum  power  is  received  when 
full  coherent  reflection  occurs,  i.e.  when  the  surface  is 
perfectly  flat  (Ohi=0).  In  such  a  condition  it  is  easy  to  verify 
that  Pnd=Pnc2  and  the  non-coherent  term  Pnc  reduces  to 
zero,  while  the  coherent  term  Pc  approaches  the  shape  of  the 
transmitted  pulse,  which  is  maximum  for  x  =  0  so  that: 

(2-2.15)  ®  i, max  (0  )  =  r,  (0  )rth 2 

which  is  a  value  consistent  with  well  know  result  for 
perfectly  flat  surfaces  reflection  coefficient.  As  the  surface 
becomes  rougher  (ah i  increases,  ahl  »  X  )  the  coherent 
component  goes  towards  zero  and  non-coherent  scattering 
becomes  dominant  (geometrical  optics  model): 


(2.2-16) 
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conventional  altimeter  mode. 

In  the  following  our  attention  will  be  devoted  only  to  the 
maximum  of  the  return  waveform. 


2.2.2  Fractal  Surface  topography  model 

Following,  recent  attempts  to  describe  the  structure  of  the 
planets  surface  by  means  of  fractals,  and  several  papers  on 
the  analysis  of  MOLA  (  Mars  Orbital  Laser  Altimeter)  data, 
from  the  Mars  Global  Surveyor,  we  tried  to  introduce  this 
model  in  our  instrument  performance  evaluation.  We  wish 
to  recall  that  MOLA  is  Laser  altimeter  whose  data  can  be 
reduced  to  topographic  height  and,  due  to  the  time  interval 


between  the  measurement,  are  spaced  of  about  300  m.  In 
particular  it  was  outlined  (s.  e.g.  Fig.2a)  from  [7])  that  Mars 
surface  point  to  point  slope  seems  to  display  a  power  law 
behavior  vs.  the  lag  distance  typical  in  fractal  profile  and  the 
parameter  that  drives  the  scaling  behavior  is  called  Hurst 
exponent  H.  Moreover  by  considering  larger  values  of  lag 
distance  we  have  obtained  a  best  fitting,  by  experimental 
analysis,  with  the  classical  stationary  Gaussian  amplitude 
model  correlated  with  exponential  correlation  function 


(2.2-17) 
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Fractal  correlation  function  [5]: 
(2.2-18) 
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where  for  stationary  functions  ks=2  [5]. 

Removing  the  average  trend  from  the  profiles  we  have 
planned  to  obtain  the  mean  value  of  the  RMS  height  from 
Mola  data 


(2.2-19)  ,=i 

The  Hurst  coefficient  can  be  estimated  trough  the  Allan 
variance  (Gallant)  (semivariogram)  that  can  be  computed 
trough  the  N  profile  samples  according  to  (Ax  =  k  -Ax0— >k) 

(2.2-20)  v(&)  ([z(x)  —z(x +Ar)]2)  =VfAv()/^-  =+(1)-^— ) 
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so  that 
(2.2-21) 
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Thence  trough  a  linear  best  fitting,  performed  in  a  proper 
lag  distance,  of  the  logarithm  of  the  Allan  variance,  H  and 
n(l)  (and  the  slope)  values  for  each  subsection  can  be 
estimated. 

Moreover  the  same  profiles  and  the  same  sections  can  be 
processed  to  obtain  the  correlation  function  according  to: 


(2.2-22) 


p{k) 


X  6(0-  z0  + k )) 

_ i= 1 _ 

1  N-k  i  N 

(z(.)y  +  |£  (z(/))* 

^  i=\  ^  i=\+k 


To  model  the  correlation  function  with  the  analytical  model 
(s.(2.2-18))  we  can  write: 

(2.2-23) 


log(ln  — 1— )  =  2H  log  k  -  2H  log  kc 
p(k) 

Again  from  a  linear  best  fitting  in  a  logarithm  scale 
estimation  of  the  parameters  H  and  kc  can  be  obtained. 

In  order  to  validate  the  selected  algorithms  we  have  selected 
a  surface  region  (3.5-4  S,  112.5-113  E)  of  10*10  Km  (due 
to  the  preliminary  information  of  the  correlation  length  and 
M ARSIS  spatial  resolution).  Taking  into  account  of  the 
behaviour  of  fig. 3a),  by  best  fitting  we  have  obtained  the 
flat  surface  of  fig.3b),  so  that  on  the  difference  (s.  fig. 3c)) 
we  have  applied  the  previous  estimation  technique  (for  each 
selected  profile  in  column  and  line  direction,  performing 
proper  mean  values)  obtaining  the  fractal  estimated 
parameters  shown  always  in  fig. 3.  We  have  obtained  H 
values  near  1,  so  that  the  geometric  optical  model  can  be 
applied. 

The  previous  approach  was  applied  to  Guesv  crater  (13-16 
S,  173-177  E,  characterized  by  big  slope  around  the  crater) 
of  170*250  Km2,  in  all  (overlapped  in  column  and  line 
direction  of  5  km)  sub  regions  of  10*10  Km2  ,  obtaining  the 
results  of  fig.4. 

We  are  arrived  to  the  conclusion  that  the  well  known 
backscattering  models  under  Kirchhoff  formulation  can  be 
applied  [3] [4]. 

Moreover  we  wish  recall  the  following  range  of  surface  size 
involved  in  MARSIS  Altimeter  by  considering  pulsed 
limited  region  and  cross  track  region  over  imposed  to 
subsurface  region  to  be  investigated: 


Radar  range  resolution  =  r 

RPL=pkRp^Rp=^^-«h^RPL=%-\5km 


Pulsed  Limited  radius 

R§  =^2SRp  S=ir — >  depth  <5  =  5km  ~^R§=44 —90km 

Cross-track  range,  from  nadir  pointing,  displacement  of  the 


surface  clutter  region  referred  to  subsurface  depth  8  : 
the  permittivity  constant  of  the  first  layer  will  be  taken  into 
account  in  8. 
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where  the  cross-track  angle  and  size  are  referred  to  the 
circular  rings. 

Concerning  the  backscattering  from  the  surface,  pulse 
limited  region,  (taking  into  account  that  we  are  in  side  or 
outside  Fresnel  region  and  the  small  observation  angle 
referred  to  nadir  pointing  )  we  wish  recall  the  cross-section 
[3] [1 1],(  related  to  the  maximum  value  of  eq.(2.2-15)): 

•  nadir  observation 
(2.2-24) 
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where  Le  (L^  )  is  the  backscattering  region  size  defined  in 
along  (cross-)  track  observation  direction. 

•  cross  track  observation,  in  the  region  related  to 
subsurface  investigation  (depth  8  ),: 

(2.2-25) 
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as  worst  case  we  can  assume  the  maximum  value  of  the 
OCcls  in  the  region  under  investigation,  in  order  to  evaluate 

the  maximum  level  of  the  clutter  involving  the 
corresponding  subsurface  region. 

Moreover  where  the  H  ^1,  taking  into  account  that,  by 
considering  small  angle  referred  to  nadir  observation,  we 
can  use  always  the  previous  eq.  where  the  slope  ms 
becomes: 


(2.2-26) 
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where  s (k)  is  the  value  of  the  rms  slope  evaluated  between 
points  separated  by  a  distance  equal  to  the  transmitted 
wavelength: 
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Moreover  we  can  notice  different  statistical  properties  in  x 
and  y  direction  in  some  regions  of  Guesv  crater,  e.g.  crater 
edge  as  shown  in  fig.4,  where  the  subsurface  investigation 
is  not  easy,  due  to  the  particular  slope  in  x  direction.  In  this 
case,  in  order  to  estimate  the  backscattering,  facets  methods 
should  be  applied. 

3  CONCLUSIONS 

In  order  to  evaluate  the  surface  backscattering,  the 
following  procedure  (by  MOLA  data)  can  be  applied: 

•  Selection  of  the  10*10  km2  regions,  overlapped  of  5 
Km  in  the  perpendicular  directions 

•  Flat  surface  estimation  for  each  region  by  best  fitting 

•  With  reference  to  the  flat  surface  for  each  region  must 
be  estimated: 

Zi  mean  value  of  the  surface  height 
axi,  ayi  surface  slopes  in  x  and  y  direction 

•  With  reference  to  the  difference  between  input  MOLA 
data  and  the  flat  surface  for  each  region  must  be 
estimated: 

Gi  rms  surface 
Li  correlation  length 
Hi  Hurst  coefficient 

The  previous  estimation  will  be  obtained  in  the  orbit  motion 
direction  (mean  value  of  multi  profiles  selected  by 
considering  the  radar  spatial  resolution)  or  the  same 
procedure  will  be  applied  a  priori  in  two  selected 
perpendicular  direction  and  finally  a  proper  or  the  mean 
value  of  the  previous  results  will  be  available. 

If  the  estimated  correlation  length  is  bigger  of  the  size  of 
selected  region  (10  km),  facets  development  (accuracy  best 
of  A/4)  should  be  applied. 
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